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Within the framework of a holographic dual model of QCD, we develop a formalism for calculating 
form factors of vector mesons. We show that the holographic bound states can be described not only 
in terms of eigenfunctions of the equation of motion, but also in terms of conjugate wave functions 
that are close analogues of quantum-mechanical bound state wave functions. We derive a generalized 
VMD representation for form factors, and find a very specific VMD pattern, in which form factors 
are essentially given by contributions due to the first two bound states in the Q 2 -channel. We 
calculate electric radius of the p-meson, finding the value (r?)c = 0.53 fm 2 . 



Introduction. - The AdS/CFT correspondence [l| 
conjectures equivalence of gravity theory on the Anti 
de Sitter space AdSs and a strongly coupled four- 
dimensional (4D) conformal field theory (CFT). The cor- 
respondence states that for every CFT operator 0{x) 
there is a corresponding bulk field <&(x, z) uniquely de- 
termined by the boundary condition (b.c.) <f>(x,z = 0) 
at the ultraviolet (UV) 4D boundary of AdS space (x 
denotes the 4D coordinates and z stands for the fifth ex- 
tra dimension). The addition of an infrared (IR) brane 
at z = zq breaks conformal invariance in the IR region, 
and allows one to have both particles and S'-matrix ele- 
ments. Due to the holographic equivalence between the 
broken CFT and the gravitational picture, the two the- 
ories have identical spectra and identical S'-matrix ele- 
ments 0. In particular, the Kaluza-Klein modes on the 
gravity side can be interpreted as bound states in the 4D 
theory. The next conjecture is that the AdS/CFT corre- 
spondence can be extended to assert that any 5D gravity 
theory on AdS 5 is holographically dual to some strongly 
coupled, large- A^ c 4D CFT (see, e.g., 0). The goal of 
holographic models of quantum chromodynamics (QCD) 
is to find such a gravity theory for which the dual theory 
is as close to QCD as possible. 

Holographic duals of QCD based on the AdS/CFT 
correspondence have been applied recently to hadronic 
physics (see, e.g., |, |, 1, @, Q, @, 0, E3, El El, El El). 
These models are able to incorporate essential proper- 
ties of QCD such as confinement and chiral symmetry 
breaking, and have demonstrated in many cases suc- 
cess in determination of static hadronic properties, such 
as resonance masses, decay constants, chiral coefficients, 
etc. Dynamic properties (form factors) have been stud- 
ied originally within the holographic approach of Ref. 

and the connection between AdS/QCD approach 
of Refs. @, Q and the usual light-cone formalism for 
hadronic form factors was proposed in [ll[ and discussed 
in [l5| • The calculation of form factors of scalar and vec- 
tor hadrons within the approach of Ref. 0] was performed 
111 Refs. HSE3, and applied to study the universality of 
the p-meson couplings to other hadrons. The expressions 



for hadronic form factors given in Refs. 0,EII,E1 have 1 
an expected form of z-integral containing the product 
of two hadronic wave functions and a function describ- 
ing the probing current. However, the hadronic functions 
used in Ref. [ll[ strongly differ from those in Refs. [l,E3l- 
The latter give meson coupling constants through their 
derivatives at z = and satisfy Neumann b.c. at the IR 
boundary z = z , while the functions used in Ref. [ll| sat- 
isfy Dirichlet b.c. at z = zq, and are proportional (after 
extraction of the overall z 2 factor) to the meson coupling 
constants /„ at the origin. In these respects they are 
analogous to the bound state wave functions in quantum 
mechanics, which makes possible their interpretation in 
terms of light-cone variables proposed in Ref. . 

The aim of this letter is to study form factors and 
wave functions of vector mesons within the framework of 
the holographic QCD model described in Refs. @, 0, @| 
(which will be referred to as H-model). To this end, we 
consider a 5D dual of the simplest Nf = 2 version of 
QCD to be a Yang-Mills theory with the SU(2) gauge 
group in the background of sliced AdS space, i.e., the 4D 
global SU(2) isotopic symmetry of Nf = 2 QCD is pro- 
moted to a 5D gauge symmetry in the bulk. Note, that 
the AdS/QCD correspondence does not refer explicitly to 
quark and gluon degrees of freedom. Rather, one deals 
with the bound states of QCD which appear as infinitely 
narrow resonances. The counterparts in the correspon- 
dence relation are the vector current J£(x) with confor- 
mal dimension A = 3 (in QCD, it may be visualized as 
q(x)t a 7 l _ l q(x) ), and the 5D gauge field A®(x, z). 

We start with recalling the basic elements of the anal- 
ysis of two-point functions (J J) given in Refs. @, 0], 
and introduce a convenient representation for the 
A-ficld bulk-to-boundary propagator V(p, z) based on the 
Knescr-Sommcrfcld formula [18| that gives V(p, z) as an 
expansion over bound state poles with the z-dcpendcncc 
of each pole contribution given by "tp wave functions" , 
that are eigenfunctions of the 5D equation of motion 
with Neumann b.c. at the IR boundary. Then we study 
the three-point function (J J J) and obtain expression for 
transition form factors that involves ip wave functions 
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and the nonnormalizablc mode factor J{Q, z). We write 
the latter as a sum over all bound states in the chan- 
nel of electromagnetic current, which gives an analogue 
of generalized vector meson dominance (VMD) repre- 
sentation for hadronic form factors. As the next step, 
we introduce "</> wave functions" that strongly resemble 
wave functions of bound states in quantum mechanics 
(they satisfy Dirichlet b.c. at z = zq, and their values 
at z = give bound state couplings g<hf n /M ni i.e., they 
have the properties necessary for the light-cone interpre- 
tation of AdS/QCD results proposed in Ref. [Hj]). We 
rewrite form factors in terms of <f> functions, formulate 
predictions for p-meson form factors, and analyze these 
predictions in the regions of small and large Q 2 . 

The p-meson electric radius is calculated, and it is also 
shown that H-model predicts a peculiar VMD pattern 
when two (rather than just one) lowest bound states in 
the Q 2 -channel play the dominant role while contribu- 
tions from higher states can be neglected. This double- 
resonance dominance is established both for the p-mcson 
form factor F(Q 2 ) given by the overlap of the -0-wave 
function (here we confirm the results obtained in Ref. [l6[ 
for the p-meson form factor considered there) and for the 
form factor T(Q 2 ) given by the overlap of the </>-wave 
function. Finally, we summarize our results. 

Two-Point Function. - Our goal is to analyze form 
factors of vector mesons within the framework of the 
holographic model of QCD based on AdS/QCD corre- 
spondence. As a 4D operator on the QCD side, we take 
the vector current J^(x) = q(x)^ ^t a q(x) , to which cor- 
responds a bulk gauge field A a M (x,z) whose boundary 
value is the source for J£(x). We follow the conventions 
of the H-model 0] , with the bulk fields in the background 
of the sliced AdSs metric 



ds 2 
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{r\ iiv dx ii dx v - dz 2 ) , < z < z , (1) 



where r?^ = Diag (1, -1, -1, -1), and z ~ I/Aqcd is 
the imposed IR scale. The 5D gauge action in AdS.5 
space, corresponding to A\j{x^ z), is 



S A ds = -^iJd 4 xdz ^gTi (F MN F MN ) 



(2) 



where F MN = d M A N -d N A M -i[AM, A N ], A M = t a A a Ml 
(t a e SU{2),a = 1,2,3) and M, N = 0, 1,2,3, z. Since 
the vector field A a M (x, z) is taken to be non-Abelian, the 
3-point function of these fields in the lowest approxima- 
tion can be extracted directly from the Lagrangian. 

Before calculating the 3-point function, we recall some 
properties of the 2-point function discussed in 0- Con- 
sider the sliced AdS space with an IR boundary at 
z = zq and UV cutoff at z = e (taken to be zero at 
the end of the calculations). In order to calculate the 
current-current correlator (or 2-point function) using the 
AdS/CFT correspondence, one should solve equations 



of motion, requiring the solution at the UV boundary 
(z = 0) to coincide with the 4D source of the vector cur- 
rent, calculate 5D action on this solution and then vary 
the action (twice) with respect to the boundary source. 
The task is simplified when the A z = gauge is im- 
posed, and the gauge field is Fourier-transformed in 4D, 
Afj,(x, z) => Apip, z). Then 



A„(p,z) = A„(p) 



V(p,e) 



(3) 



where A^(p) is the Fourier-transformed current source, 
and the 5D gauge field V(p, z) is the so-called bulk-to- 
boundary propagator obeying 



zd z [-d z V{p,z)) +p 2 V(p,z) = 



(4) 



The UV b.c. A fJi (p,e) — A^(p) is satisfied by construc- 
tion. At the IR boundary (when z = zo), we follow 
Ref. [?| (see also Ref. [HI) and choose the Neumann b.c. 
d z V(p, zq) = which corresponds to the gauge invariant 
condition i 7 ' A12 (x, zo) = 0. Evaluating the bilinear term of 
the action on this solution leaves only the UV surface 
term 



? (2) 
'AdS 



^2 / ^7^(p)A,(p) 



WSJ (27T) 



ld z V(p, z) 
z V(p,e) 



(5) 



The 2-point function of vector currents is defined by 



d A x e^{J«(x)4(0)) = 5 ab n^(p)Z(p 2 ) , (6) 



where II^p) = (?7 M „ — p^p^/p 2 ) is the transverse pro- 
jector. Varying the action ([5]) with respect to the bound- 
ary source produces 



5V) 



1 fl d z V(p,z) 
z V(p,e) 



9% 



(7) 



(To get the tensor structure of ((6|) by a "naive" variation, 
one should change A^A^ -> A»U U y(p)A v in Eq. ©). 

It is well known (see, e.g., [3, [Ta]) that two linearly 
independent solutions of Eq. ^ are given by the Bcsscl 
functions zJi(Pz) and zYi(Pz), where P = \fp 2 . Taking 
Neumann b.c. for V(p, z), one obtains 



V(p, z) = Pz 
and, hence, 



Y (Pz )J 1 (Pz) - J (Pz )Y 1 (Pz) 



S(P 2 ) 



np 

2g? 



2 r 



5 L 



Y (Pz) - J (Pz) 



Y (Pz ) 



This expression is singular as e 
1 

2ff 5 



■h(Pz ) 
0: 



S(p 2 



2 P 2 ln(p 2 e 2 ) + 



(8) 



(9) 



(10) 
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By matching to QCD result for J° = <?7 M i a <? currents one 
finds flf = 12ir 2 /N c (cf. d). 

The two-point function S(p 2 ) has poles when the de- 
nominator function Jq(Pzq) has zeros, i.e., when Pzq co- 
incides with one of the roots jo >n of the Bessel function 
Jo(x). These poles can be explicitly displayed by incor- 
porating the Kneser-Sommerfeld expansion (l8j 

Y (Pz ) J (Pz) - J (Pz Q )Y (Pz) 



MPzo) 



Joho,nZ/z ) 



* ^[Hl»,n)?{P 2 4-lln) 



(11) 



valid for z < zq (the case we are interested in). Taking 
formally z = gives a logarithmically divergent series 
reflecting the In e singularity of the z = e expression. 
Thus, some kind of regularization for this divergency of 
the sum is implied. Under this assumption, 



Y( ^ ^ V- [ J l(70.n)]" 2 



(12) 



where M n = 7o,n/zo- Hence, the 2-point correlator of the 
H-model has poles when P coincides with one of M n 's. 
Given that the residues of all these poles are positive, the 
poles may be interpreted as bound states with M n 's being 
their masses. The coupling f 2 with which a particular 
resonance contributes to the total sum is determined by 

f n = lim {(p 2 -M 2 )Z(p 2 )} . (13) 

P ->M* 



(0|^l/°fi) = $ ab fn e n for the vector meson decay con- 



This prescription agrees with the usual definition 
(0\J a M) = S ab f n e 
stants. In our case, 

f 

J n 



2Ml 



(14) 



Three-Point Function. - Consider now the trilincar 
term of the action calculated on the V(q, z) solution: 



,5' 



(3) 
AdS 



^abc 

w 5 



z ° d? 

d 4 x I — (d^)A^»A v ' c . (15) 



A naive variation gives the result for the 3-point cor- 
relator (Ja{Pi)Jb{—P2)Jc{q)) that contains the isotopic 
Levi-Civita tensor e a bc, the dynamical factor 

W(p 1 ,p 2 ,q)= / —-777 TTTT r r, (16) 



; z V(pt,e) V(p2,e) V(q,e) 
and the tensor structure 
T a ^ = rj^(q - Pl f- V ^(P2 + q) a + V^iPi + V*T 

familiar from the QCD 3-gluon vertex amplitude. Restor- 
ing the transverse projectors IT"" (pi), etc. one can con- 
vert it into 

T q/3ai = ij a/ j(pi + p 2 ) M + - ri 0fl q a ) . (17) 



For the factors corresponding to the hadronized channels, 
the Kncscr-Sommcrfcld expansion (jll[) gives 



V(p,z) 



V(p,z) 



V(p,Q) 

where p equals p\ or pi , and 

V2 



P 2 -M 2 



i>n(z) 



zJ t {M n z) 



(18) 



(19) 



is the "-0 wave function" obeying the same equation of 
motion ([4| as V(p, z) (with p 2 = M 2 ), satisfying the b.c. 



^„(o) = o , dMz ) = o 

and normalized according to 

" z ° dz 



Wn{z)Y 



1 . 



(20) 



(21) 



One remark is in order here. Since the u ip wave func- 
tions" vanish at the origin and satisfy Neumann b.c. at 
the IR boundary, it is impossible to establish a direct 
analogy between ^ n (z)'s and the bound state wave func- 
tions in quantum mechanics. For the latter, one would 
expect that they vanish at the confinement radius, while 
their values at the origin are proportional to the coupling 
constants /„. 

Taking a spacelike momentum transfer, q 2 = —Q 2 for 
the V/V factor of the EM current channel gives 



J(Q,z) = Qz 



K 1 {Qz) + I 1 (Qz) 



K (Qz ) 



Io{Qzo) 



(22) 



the non-normalizable mode with Neumann b.c. (see also 
Ref. [3]). This factor can also be written as a sum of 
monopole contributions from the infinite tower of vector 
mesons: 



Q 2 + Ml 



(23) 



This decomposition, discussed in Ref. [16j, directly fol- 
lows from Eq. (|18p . Incorporating the representation for 
the bulk-boundary propagators given above we obtain 



t(pIpIq 2 )= £ 



fnfkFnk(Q 2 ) 



Ml) 



(24) 



where T(p 2 ,plQ 2 ) 

Fnk(Q 2 ) = 



c=i (ri-M^(p 

W(pi,p2,q)/gl, and 
z ° dz 

— J(Q,z)Mz)Mz) (25) 



correspond to form factors for n —> k transitions. This 
expression was also written in Ref. (lij for form factors 
considered there. 
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Wave functions. - The formulas obtained above us- 
ing explicit properties of the Besscl functions in the form 
of Kneser-Sommerfeld expansions, can also be derived 
from the general formalism of Green's functions. In par- 
ticular, the Green's function for Eq. can be written 
as 



1pn{z)ll) n (z') 



(26) 



where ip n {zYs are the normalized wave functions (|19j) 
that satisfy the Sturm-Liouvillc equation (H} with 
p 2 = M 2 and Neumann b.c. (|20[) . As discussed in 
Ref. [||, the bulk-to-boundary propagator is related to 
the Green's function by 



V(p,z') = - 



d z G(p; z, z') 



(27) 



and the two-point function T,(P 2 ) is obtained from the 
Green's function using Eqs. f7|. (|27l) 



£(P 2 ) = -j 

95 



-d z G(p;z,z') 



(28) 



Accordingly, the coupling constants are related to the tp 
wave functions by 



fn 



1 



95 iz 



dzlpn(z) 



(29) 



2=0 



(cf. [6j, UM). In view of this relation, it makes sense to 
introduce u d> wave functions" 



4>n{z) 



1 



M n z 



9 z 1pn(z) = 



V2 



ZoJl(lO,n) 



J (M n z) , (30) 



which give the couplings g5f n /M n as their values at the 
origin. In this respect, the "(f) wave functions" are analo- 
gous to the bound state wave functions in quantum me- 
chanics. Moreover, these functions satisfy Dirichlct b. c. 
i'nizo) = and are normalized by 



dzz\(j> n {z)\' 



1 



(31) 



which strengthens this analogy. However, the elastic form 
factors F nn (Q 2 ) are given by the integrals 



Fnn{Q ) 



dz 



J(Q,z)\Mz)f 



(32) 



involving ip rather than </> wave functions. In fact, due 
to the basic equation (|1]), tp n (z) wave functions can be 
expressed in terms of 4> n (z) as 



ipn(z) 



M r , 



d z (/)n(z) , 



(33) 



and we can rewrite the form factor integral as 



dzzJ(Q,z)\ct) n (z)\' 



(34) 



1 f*° 

+ TT I dz(/) n (z)ip n (z)d z J(Q,z) 

Win. .In 



Note, that the nonnormalizablc mode 
1 



■ d z J(Q, z) = -Q' 



K (Qz) - I (Qz) 



Ko(Qz ) 



h{Qzo) 



(35) 



corresponds to equation whose solutions are the functions 
Jo{M n z) satisfying Dirichlct b.c. at z = zq. Expressing 
4> n (z) in terms of d z ip n (z), integrating |f/j„(z)| 2 by parts 
and using equation (01 for J{Q, z) gives 



Fnn 



dzzJ(Q,z) |0„(z)p 



(36) 



Q 2 
2M 2 



— J(Q,z)\Mz)f 



The second term contains the original integral for 
Fnn(Q 2 ), and we obtain 



Fnn(Q ) 



1 



1 + Q 2 /2M 2 



dzzJ{Q,z)\4> n {zf 



(37) 

Notice, that the normalizable modes 4> n (z) in this expres- 
sion correspond to Dirichlct b.c, while the nonnormal- 
izable mode ^T(Q, z) was obtained using the Neumann 
ones. 

Thus, we managed to get the expression for F nn {Q 2 ) 
form factors that contains <fi instead of tp wave functions. 
However, it contains an extra factor 1/(1 + Q 2 /2M 2 ), 
which brings us to the issue of different form factors of 
the p-mcson and kinematic factors associated with them. 

Form Factors. - Our result (|25|) contains only one 
function for each n —* k transition, in particular F nn (Q 2 ) 
in the diagonal case. However, the general expression 
for the EM vertex of a spin-1 particle of mass M can be 
written (assuming P- and T-invariance) in terms of three 
form factors (see, e.g., 0, our G 2 is theirs G% — G\): 



( P + (p2,e')\j£ M (Q)\P + (Pi,t)) 
= -e' e a [r ] ^(p^+p^G 1 (Q 2 ) 

+(tt"V " ^V)(Gi(Q 2 ) + G 2 (Q 2 )) 



(38) 



1 

M 



Comparing the tensor structure of this expression with 
(fTT)) . we conclude that H-model predicts (Q 2 ) = 
G { 2 ] {Q 2 ) = F nn (Q 2 ), and G { 3 n) (Q 2 ) = for form fac- 
tors G[ n) (Q 2 ) of n th bound state. It was argued (sec 
[l6|) that this is a general feature of AdS/QCD models 
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for the p-meson form factors. Since J(Q = 0, z) = 1, the 
diagonal form factors F nn (Q 2 ) in the H- model are nor- 
malized to unity, while the nondiagonal ones vanish for 
Q 2 = (the functions tp n (z) are orthonormal on [0, Zq])- 
The form factors Gi are related to electric Gc, mag- 
netic Gm and quadrupole Gq form factors by 



Gc = G\ 



G c 



1 + 



Q 2 

6M 2 
AM 2 



Gq , Gm — Gi + G2 



G3 — G9 



(39) 



For these form factors, H- model gives 
G { q\Q 2 ) = -F nn {Q 2 ) , G$(Q 2 ) = 2F nn (Q 2 ) , (40) 



g ( S\q 2 ) = [i 



6M : 



Fnn (Q ) 



For Q 2 — 0, it correctly reproduces the unit electric 
charge of the meson, and "predicts" p = Gm (0) = 2 for 
the magnetic moment and D = Gq (0)/A/ 2 = — 1/M 2 for 
the quadrupole moment, which are just the canonical val- 
ues for a pointlike vector particle [201 ] . 

Another interesting combination of form factors 



F{Q 2 ) = G 1 (Q 2 



2M- 



G 2 (Q 2 )- 



2AF 



G 3 (Q 2 



(41) 

appears if one takes the "+++" component of the 3-point 
correlator (obtained, e.g., by convoluting it with n^n^n^, 
where n 2 = 0, (npi) = l,(nq) = [la]). The H-modcl 
result ([57)) for F(Q 2 ) is particularly simple: 



•F nn(Q ) 



dzzJ{Q,z) \4> n (zf 



(42) 



Thus, it is the form factors F nn (Q 2 ) that are the most 
direct analogues of diagonal bound state form factors in 
quantum mechanics. 

Low-Q 2 behavior. - Our expression for F nn (Q 2 ) is 
close to that proposed for a generic meson form factor 
in the holographic model of Ref. ll|. There, the au- 
thors used JC(Qz) = QzKi(Qz) as the ^-channel factor. 
Indeed, the difference between J{Q, z) and K,(Qz) is ex- 
ponentially small when Qzq 3> 1, but the two functions 
radically differ in the region of small Q 2 , where the func- 
tion IC(Qz) displays the logarithmic branch singularity 



K{Qz) 



z 2 Q 2 



l-2 7B -rn(QV/4) 



-0(Q 4 ) , (43) 



that leads to incorrect infinite slope at Q 2 = 0. To im- 
plant the AdS/QCD information about the hadron spec- 
trum in the g-channel one should use J{Q, z) that corre- 
sponds to a tower of bound states in the g-channel. The 
lowest singularity in this case is located at Q 2 = —M\. 



Since it is separated by a finite gap from zero, the form 
factor slopes at Q 2 = are finite. 

To analyze the form factor behavior in the Qzq <C 1 
limit, we expand 



z 2 Q 2 



J(Q,*)\c 



'zo<«l 



= 1 - 



1 - In- 



0{Q A ) . (44) 



As expected, the result is analytic in Q 2 . For the low- 
est transition (i.e., for the p-meson form factor), explicit 
numbers are as follows: 



Fn{Q 2 ) w 1-0.692 



Q' 2 



0.633 



Q 4 



0(Q 6 ), (45) 



M 2 ' M 4 
where M = Mi = m p . Another small-Q 2 expansion 

F n (g 2 ) » 1 - 1.192 i£ + 1.229 ^ + 0(Q 6 ) , (46) 

can be either calculated from the original expression 
([32]) involving ^-functions or by dividing T\\ (Q 2 ) by 
(1 + Q 2 /2M 2 ). The latter approach easily explains the 
difference in slopes of these two form factors at Q 2 = 0. 
Finally, for the electric form factor, we obtain 



G[]\Q 2 ) m 1 - 1.359 



1.428 
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0(Q 6 ) . (47) 



M 2 M 4 
For the electric radius of the p-mcson this gives 

(r 2 ) c = 0.53 fm 2 , (48) 

the value that is very close to the recent result (0.54 fm 2 ) 
obtained within the Dyson-Schwinger equations (DSE) 
approach [2l[. Lattice gauge calculations [22j indicate a 
similar value in the m 2 — > limit. 

Vector meson dominance patterns. - Numeri- 
cally, the result 1.359/M 2 for the slope of G^ (Q 2 ) is 
larger than the simple VMD expectation 1/M 2 . In fact, a 
part of this larger value is due to the factor (1 — Q 2 /6A/ 2 ) 
relating Gq\q 2 ) and Fn{Q 2 ), which is kinematic to 
some extent. The Fii(Q 2 ) form factor, however, can be 
written in the generalized VMD representation (cf. 



Fn(Q 2 )=J2 Y 



F„ 



(49) 



with the coefficients F m ,n given by the overlap integrals 

Jl(70,m^) Ji(70,l0 



m,ll 



= 4 



^l( 7 0,m)Ji 2 (70,l) 



(50) 



apparently having a purely dynamical origin. The coef- 
ficients F TOi n satisfy the sum rule 



TCI — 1 



(51) 
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that provides correct normalization Fn(Q 2 = 0) = 1. 
Numerically, the unity value of the form factor Fn(Q 2 ) 
for Q 2 = is dominated by the first bound state that 
gives 1.237. The second bound state makes a sizable cor- 
rection by —0.239, while adding a small 0.002 contribu- 
tion from the third bound state fine-tunes 1 beyond the 
10~ 3 accuracy. Contributions from higher bound states 
to the form factor normalization are negligible at this 
precision. 

The slope of Fn(Q 2 ) at Q 2 = is given by the sum of 
Fmsi/M^ coefficients. Now, the dominance of the first 
bound state is even more pronounced: the Q 2 coefficient 
1.192/M 2 in Eq. (j46]) is basically contributed by the first 
bound state that gives 1.237/Af 2 , with small -0.045/M 2 
correction from the second bound state. Other reso- 
nances are not visible at the three-digit precision. 

Thus, for small Q 2 , H-model predicts a rather pecu- 
liar pattern of VMD for Fn(Q 2 ) (observed originally in 
Ref. [l6| for a form factor considered there): strong dom- 
inance of the first g-channcl bound state, whose coupling 
exceeds 1, with the second resonance (having the 
negative coupling i<2,ii) compensating this excess. 

Similarly, the Tn(Q 2 ) form factor has the generalized 
VMD representation with coefficients T m ,ii given by the 
overlap integrals 



J 11 



Ji(7o,m)</ 2 (7o,i) 



(52) 



Now, ^i,!! w 0.619, ^2,11 ~ 0.391, ^3.11 ~ -0.012, 
•^4,11 ~ 0.002, etc. In this case also, the value of the 
(Q 2 ) form factor for Q 2 = is dominated by the 
first two bound states. For the slope of the form fac- 
tor at Q 2 = 0, the dominance of the first bound state is 
again more pronounced: the Q 2 coefficient 0.692/M 2 in 
Eq. (|4"5")) is basically contributed by the first bound state 
that gives 0.619/Af 2 , with a small 0.074/M 2 correction 
from the second bound state and a tiny — 0.001/M 2 con- 
tribution from the third one. 

Thus, for T\i(Q 2 ), H-model gives again a two- 
resonance dominance pattern, with the coupling JF2.11 of 
the second resonance being now just somewhat smaller 
than the coupling T\_\\ of the first resonance, both being 
positive. The relation between the two VMD patterns 
follows from Eq. (j3T)) : 



1 - M 2 J2M 2 



(53) 



In particular, it gives Fi n = ^^1,11, and negative sign 
for F 2 .\i- It also determines that if higher coefficients 
•Fm.ii are small then -Fail's are even smaller. 

Large-Q 2 behavior. - Eq. (|37|) tells us that asymp- 
totically Fn(Q 2 ) is suppressed by a power of 1/Q 2 com- 
pared to -Fii(Q 2 ), which is known to behave like 1/Q 2 
for large Q 2 [VgUM- The absence of 1/Q 2 ter m in the 



asymptotic expansion for Fn(Q 2 ) means that the coef- 
ficients F m , 11 defined in Eq. (149)) satisfy the "supercon- 
vergence" relation 
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(54) 



reflecting a "conspiracy" [l(| between the poles. Writ- 
ing M m F m ^i = A m M 2 , we obtain that A 1 m 1.237, 
A 2 « -1.261, A 3 w 0.027 (our results for the ratios 
A2/A1, A3/A1 agree with the calculation of Ref. [l~6j]V 
Again, the sum rule is practically saturated by the first 
two bound states, which give contributions that arc close 
in magnitude but opposite in sign. 

In case of J-(Q 2 ), the two lowest bound states both 
give positive 0(1/Q 2 ) contributions at large Q 2 . In 
Ref. [15| . it was shown that the asymptotic normaliza- 
tion of Tn{Q 2 ) exceeds the VMD expectation Mf/Q 2 by 
a factor of 2.566. We can infer this normalization from 
the values of the coefficients F m ,ii defined in Eq. |52|) . 
Writing M^Tm, 11 = AmM?, we obtain that Ai « 0.619, 
A 2 « 2.061, A 3 w -0.150, A4 w 0.054. Note, that 
the total result is dominated by the second bound state, 
which is responsible for about 80% of the value. The 
lowest bound state contributes only about 25%, while 
the higher states give just small corrections. 

It is worth noting that the large-Q 2 behavior of 
J~n (Q 2 ) is determined by the large-Qzo form of J{Q, z): 
it can be (and was) calculated using JC(Qz), the free-field 
version of J{Q, z). As a result, the value of the asymp- 
totic coefficient (2.566 in case of ^ r n(Q 2 )) is settled by 
the sum rule 



I0i(O)| 2 



dxX 2 Ki( X ) 



(55) 

that should be satisfied by any set of coefficients T m ^\\. 
A particular distribution of "2.566" among the bound 
states is governed by the specific g-channel dynamics (in 
our case, by the choice of the Neumann b.c. for J{Q, z) 
at z = zq). Thus, in the dynamics described by J{Q, z), 
the large value of the asymptotic coefficient is explained 
by large contribution due to the second bound state. 

It was shown in Ref. [l5| that the asymptotic 1/Q 2 be- 
havior for Fn(Q 2 ) is established only for Q 2 ~ lOGcV 2 , 
and one may question the applicability of the H-model for 
such large Q 2 . The discussion of this problem, however, 
is beyond the scope of the present paper. 

Summary. - In this letter, we described the formal- 
ism that allows to study form factors of vector mesons in 
the holographic QCD model of Refs. 0,0,1 (H-model). 
An essential ingredient of our approach is a systematic 
use of the Kneser-Sommerfeld representation that explic- 
itly displays the poles of two- and three-point functions 
and describes the structure of the corresponding bound 
states by cigenfunctions of the 5D equation of motion, 
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the u i/j wave functions" . These functions vanish at z = 
and satisfy Neumann b.c. at z = zq, which prevents a 
direct analogy with bound state wave functions in quan- 
tum mechanics. To this end, we introduced an alter- 
native description in terms of "</> wave functions" that 
satisfy Dirichlet b.c. at z = zq and have finite values at 
z = which determine bound state couplings g5f n /M n . 
Thus, the 4> wave functions have the properties necessary 
for the light-cone interpretation proposed in Ref. [lH and 
discussed also in Ref. [la]. 

Analyzing the three-point function, we derived expres- 
sions for bound state form factors both in terms of ip 
and 4> wave functions, and obtained specific predictions 
for form factor behavior at small and large values of the 
invariant momentum transfer Q 2 . In particular, we cal- 
culated the electric radius of the p meson, and obtained 
the value {r 2 )c = 0.53 fm 2 that practically coincides with 
the recent result |2l| obtained within the DSE approach. 
Our result is also consistent with the m 2 — ► extrapola- 
tion of the recent lattice gauge calculation [jjj ]. 

We derived a generalized VMD representation both for 
the Fn(Q 2 ) form factor (the expression for which coin- 
cides with a model p-meson form factor considered in 
Ref. [lij]) and for the T\\{Q 2 ) form factor introduced in 
the present paper, and demonstrated that H-model pre- 
dicts a very specific VMD pattern, in which these form 
factors are essentially given by contributions due to the 
first two bound states in the Q 2 -channel, with the higher 
bound states playing a negligible role. We showed that, 
while the form factor slopes at Q 2 — in this picture are 
dominated by the first bound state, the second bound 
state plays a crucial role in the large-Q 2 asymptotic limit. 
In particular, it provides the bulk part of the negative 
contribution necessary to cancel the naive VMD l/Q 2 
asymptotics for the Fn(Q 2 ) form factor (corresponding 
to the overlap integral involving the ip functions), and it 
dominates the asymptotic l/Q 2 behavior of the J-(Q 2 ) 
form factor (given by the overlap of the <f> functions). 

A possible future application of our approach is the 
analysis of bound state form factors in the model of 
Ref. [Tit ] in which the hard-wall boundary conditions at 
the z = zq IR boundary are substituted by an oscillator- 
type potential. This model provides the M 2 ~ riA 2 
asymptotic behavior of the spectrum of highly excited 
mesons, which is more consistent with the semiclassical 
limit of QCD jH] than the M 2 ~ n 2 A 2 result of the 
H-model. 
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